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ABSTRACT 

The goal of the Gemini Deep Deep Survey (GDDS) is to study an unbiased sample of K < 20.6 galaxies in the 
redshift range 0.8 < z < 2.0. Here we determine the statistical properties of the heavy element enrichment in 
the interstellar medium (ISM) of a subsample of 1 3 galaxies with 1 .34 < z < 1 .97 and UV absolute magnitude 
^2000 < -19.65. The sample contains 38% of the total number of identified galaxies in the first two fields 
of the survey with z > 1.3. The selected objects have colors typical of irregular and Sbc galaxies. Strong 
[O II] emission indicates high star formation activity in the H II regions (SFR ~ 13- 106 M yr _1 ). The high 
S/N composite spectrum shows strong ISM Mg II and Fe II absorption, together with weak Mn II and Mg I 
lines. The Fe II column density, derived using the curve of growth analysis, is log/VFeii = 15.54^(j |3. This is 
considerably larger than typical values found in damped Lya systems (DLAs) along QSO sight lines, where 
only 10 out of 87 (~ 11%) have log/Vpeii > 15.2. High Fe II column densities are observed in the z = 2.72 
Lyman break galaxy cB58 (log/Vpen ~ 15.25) and in gamma-ray burst host galaxies (log/Vpe!! ~ 14.8- 15.9). 
Given our measured Fe II column density and assuming a moderate iron dust depletion (<5p e ~ 1 dex), we derive 
an optical dust extinction Ay ~ 0.6. If the H I column density is log/Vni < 21.7 (as in 98% of DLAs), then 
the mean metallicity is Z/Z Q > 0.2. The high completeness of the GDDS sample implies that these results are 
typical of star-forming galaxies in the 1 < z < 2 redshift range, an epoch which has heretofore been particularly 
challenging for observational programs. 

Subject headings: cosmology: observations - galaxies: abundances - galaxies: ISM 



1. INTRODUCTION 

The interstellar medium (ISM) is an important component 
in galaxies, not only in terms of mass, but also because it car- 
ries information on past generations of stars, and provides the 
main fuel for the new generations of stars. In the local Uni- 
verse the cold ISM (T < 1000 K) is -15-40% of the total 
mass in stars (Fukugita et al. 1998; Zwaan et al. 2003), and 
at higher redshifts it is a major constituent of the Universe 
(~ 2 times higher than in the local Universe at 1.5 < z < 2.0; 
Ston-ie-Lombardi & Wolfe 2000; Peroux et al. 2003). In the 
cold ISM, heavy elements - the footprint of the stellar activity 
- can most easily be investigated through the analysis of UV 
absorption lines. 

Apart from very detailed analysis in Galactic ISM (Savage 
& Sembach 1996), in the local Universe exploring the UV 
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absorption lines become possible only for a few bright star- 
forming galaxies (Heckman et al. 2001; Aloisi et al. 2003) 
due to limited capabilities of UV satellites. At higher red- 
shifts, bright QSO sources, whose sight lines cross interven- 
ing galaxies, allow the detection of absorption lines with UV 
satellites and (at z > 2.5) ground-based telescopes. QSO 
damped Lyman-a systems (QSO-DLAs) have probed the 
ISM in more than 100 galaxies at 0.0 < z < 4.5 (Prochaska et 
al. 2003a). These data may be biased, however, to galaxies or 
regions of galaxies that are dust- and metal-poor since dusty, 
metal rich systems are likely to obscure or, at least, extinct the 
background QSO. 

Indeed, we have indications that the ISM in other high 
redshift objects is chemically more evolved. For instance, 
gamma-ray burst (GRB) afterglows are associated with re- 
gions of intense star formation (Bloom, Kulkarni, & Djorgov- 
ski 2002; Djorgovski et al. 2003). They are - for a short 
time - much brighter than QSOs (~ 1000 times brighter in 
X-ray fluxes) and can be detected even if embedded in a very 
dense medium. The analysis of UV absorption lines in 4 GRB 
spectra (GRB-DLAs) shows high column density of metals 
in the ISM of the host galaxies (Savaglio, Fall, & Fiore 2003; 
Vreeswijk et al. 2003). In addition, large surveys of z > 2.5 
Lyman-break galaxies (LBGs, Steidel et al. 2003) show pre- 
dominantly high equivalent widths (EWs) of ISM metal ab- 
sorption lines (Pettini et al. 2002; Shapley et al. 2003). 

The direct study of cold ISM in galaxies requires reason- 
able S/N of the rest-frame UV continuum, and this is gener- 
ally not sufficiently strong for our observational capabilities, 
even when redshifted to the optical at high z. Additionally, 
systematic problems arise from imperfect sky subtraction and 
CCD fringe removal. As a consequence, surveys have mainly 
targeted z > 2.5 galaxies where the Lyman-break technique 
starts to be effective, or strong emission line objects, observ- 
able in the optical up to z — 1. The 1 < z < 2.5 range is often 
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referred to as the "redshift desert" for lack of discoveries. This 
interval is also very problematic for studies of galaxy evolu- 
tion, since it apparently spans a major epoch of galaxy build- 
ing: 50 - 75% of the stellar mass in galaxies was in place 
by z ~ 1, whereas at 2.5 < z < 3.5 the Lyman break galax- 
ies contain only ~ 3 - 20% of the present day stellar mass 
(Dickinson et al. 2003; Fontana et al. 2003). Recently, the 
K20 survey (Cimatti et al. 2002a) with 82 spectroscopically 
confirmed K < 20 galaxies at 1 < z < 2 have started to un- 
veil some aspects of galaxy evolution in the redshift desert 
(Cimatti et al. 2002b; Pozzetti et al. 2003). However, the 
number of galaxies drops to 17 if the 1.3 < z < 2 interval is 
considered. 

To improve optical spectroscopy, a new observing mode 
known as Nod & Shuffle (N&S), has been recently devel- 
oped (Glazebrook & Bland-Hawthorn 2001; Cuillandre et 
al. 1994). During N&S observations, the telescope is rapidly 
switched between object and sky positions ("nodding"), while 
the charge is "shuffled" between different parts of the CCD. 
Because both the sky and objects are observed through the 
same optical path, slits, and pixels, the effects of fringing are 
eliminated and the sky subtraction greatly improved. As a 
result, N&S allows very long integration times before sys- 
tematic effects start to dominate, ^10 times longer than is 
practical with conventional spectroscopy. 

We have begun a project to observe galaxies in the redshift 
desert (the Gemini Deep Deep Survey, GDDS, Abraham et 
al. in preparation) that uses the N&S technique on the multi- 
object spectrograph at the Gemini 8m telescope. The ~ 10 5 
seconds exposures allow us to reach very faint magnitudes 
(limiting magnitude Iab = 24.8 and K = 20.6, more than half a 
magnitude deeper than K2Q). Both raw and fully reduced data 
from the GDDS will be publicly available 14 . In this paper we 
present first results on the heavy element enrichment in the 
cold ISM of a subset of galaxies at 1.3 < z < 2.0, based on 
observations of the first two fields. We estimate the column 
densities of Fe II, Mn II, and Mg II using the curve of growth 
analysis, and compare our results with what is already known 
of the ISM over a larger redshift scale. 

The paper is structured as follows: in §2 we present the 
data set, the sample selection is described in §3, in §4 we 
report [O II] line flux measurements, in §5 we describe the 
method used to combine the individual spectra, in §6 we mea- 
sure column densities of metals, in §7 we report results of 
other surveys, in §8 we discuss the general properties of the 
selected sample, and in §9 we summarize the main conclu- 
sions. Throughout the paper we adopt a h = H / 100 = 0.7, 
Q.M = 0.3, Ha = 0.7 cosmology. 

2. GDDS DATA 

The spectra presented in this paper were obtained during 
three observing runs at the Gemini North telescope in the pe- 
riod 2002 August-October. The two fields SA22 (a2000 = 
22 A 17 m 20", <5 2 ooo = +00°21'30") and NDWS 15 (a 2 ooo = 
02 /l 10 m 00 1 , (5 20 oo = -04°30'00") were selected from the Las 
Campanas Infrared (LCIR) Survey (McCarthy et al. 2001; 
Chen et al. 2002; Firth et al. 2002; Chen et al. 2003), and 
observed with the Gemini multi-object spectrograph (GMOS, 
5.5' x 5.5' field of view; Murowinski et al. in preparation; 
Hook et al. 2003). The "queue scheduled" mode at the Gem- 

14 Data can be downloaded at the GDDS homepage 
http ://w w w. ociw.edu/lcirs/gdds .html 

13 NDWS: NOAO Deep Wide-Field Survey. 



ini observatory ensured that data were taken only when the 
weather conditions were as requested. As a result, during our 
observations the seeing was ^ 0.85" in / band for most of the 
time. The LCIR Survey photometry (VIzHK filters for field 
SA22 and BVRIzK filters for field NDWS) was used to prese- 
lect z > 0.8 objects in the two fields with the photometric red- 
shift technique. Two GMOS masks were used for field SA22 
(exposure times 48600 and 90000 seconds) with 111 slits to- 
tal, and one mask with 60 slits for field NDWS (exposure time 
75600 seconds). The total number of objects observed was 
then 171. The observations were taken using 0.75" x 2.2" 
and 0.75" x 2" slit apertures. The comparison with similarly 
faint IR galaxies at z > 1 (Yan et al. 1998) suggests that the 
slit apertures allow to collect more than 50% of the galaxy 
flux. The 2D data reduction was performed using IRAF 16 and 
Gemini/GMOS tools 17 . For the ID spectrum extractions we 
used iGDDS, a package specifically written with the toolkit 
Cocoa under Mac OS X (details in Abraham et al. in prepara- 
tion). Uncertainty in the wavelength calibration is of the order 
of 1.5 A. The wavelength calibrated spectra were corrected 
for the telluric absorption and flux calibrated. Most of the 
observations were carried out near transit and no corrections 
for differential atmospheric refraction (which is small at these 
wavelengths) were applied. To confirm that the N&S tech- 
nique was working as expected, sky residuals were measured 
from some strong lines in a stack of 60 spectra. The residuals 
are of the order of 0.05 - 0. 1 %, as predicted by Glazebrook & 
Bland-Hawthorn (2001). The final spectral coverage in each 
spectrum is within the 5000- 10000 A interval, depending 
on the position of the slit in the mask, and the dispersion is 
AA ~ 3.5 A per pixel. The resolution element, determined by 
the slit width, is 4.5 pixels or FWHM ~ 16 A, equivalent to 
960-460 km s , from the blue to the red end of the spectral 
range. The resolution was also measured from unresolved sky 
emission lines to be FWHM ~ 17- 18 A, slightly larger than 
the nominal resolution. This is not unexpected since the sky 
emission is extended over the whole slit aperture. 

The spectra were extracted both using the linear and the op- 
timal extraction method. In parallel the noise spectrum was 
created for each object using the object and sky extracted 
signal, and the read-out-noise of the instrument (3-3.5 e~, 
depending on the CCD). Redshifts were found from simple 
eye inspection of line positions. Manual redshifting (possible 
thanks to the small number of spectra) was done because of 
the ability of the eye to reject extraneous data artifact features 
in the 2D spectrum. Among the 171 objects observed, 7 have 
been identified as stars, 4 as AGNs, 18 have tentative redshift, 
18 have no redshift. The remaining 124 extragalactic objects 
have secure redshift. 

3. SAMPLE SELECTION 

To measure metal enrichment in galaxies in the 1 < z < 
2.5 range, we use absorption lines in the rest-frame interval 
2200-2900 A. This is a particularly interesting region be- 
cause it contains 7 Fe II transitions, with very different os- 
cillator strengths, together with the strong Mg II doublet and 
the Mg I line. The Mg II doublet is normally very strong and 

16 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 

17 For more information visit 

http://www.gemini.edu/sciops/instruments/gmos/gmosIndex.html 
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Table 1 
Selected galaxy sample 
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"Redshift of the ISM absorption lines. 
b Redshift of the [O II] emission line. 
'"Velocity difference between Zroil] andzisM- 

d AB absolute magnitude in a synthetic box filter centered at A = 2000 A. The accuracy is ~ 0.2 magnitudes (see text). 
c The [O II] fluxes have been measured from the spectra, and no aperture correction has been applied. Errors include only photon 
counting statistics. 

f SFRs are estimated assuming SFR[ n](M a yr -1 ) ~ 2.54 x t[on] x 10 -4 ' erg s cm -2 . 
£SFRs are estimated assuming SFR 2 ooo(M s yr 1 ) ~ 1 .25 x L 2 ooo x 10~ 28 erg s -1 Hz~' . 



-16 



-1! 



-20 - 



-22 




1.8 2 

redshift 



FIG. 1 . — AB M2000 absolute magnitude of GDDS 1.3 < z < 2.4 galaxies 
as a function of redshift. Fluxes are extracted from a 2 arcsec aperture. Open 
and filled dots are galaxies with secure spectroscopic redshift, respectively. 
Filled dots are the selected galaxies with M2000 < —19.65. Stars are objects 
with tentative redshift. Crosses have no spectroscopic redshift; in these cases 
we used the photometric redshift. The selected sample contains between 25 
and 38% of the total number of galaxies identified in GDDS with 1.3 < z < 
2.4 (depending whether galaxies with tentative or photometric-redshift only 
are included in the count or not). 



highly saturated but it is very useful in determining galaxy 
redshifts. The useful wavelength range for most GDDS spec- 
tra is limited to A = 5000-9800 A, or the rest-frame interval 
A = 2200-2900 A for galaxies at 1.3 < z < 2.4. In the total 
sample of 124 objects with spectroscopic redshift in the two 
fields, 34 galaxies are in this z range. The number of objects 
with tentative or photometric-only redshifts in 1 . 3 < z < 2 .4 is 
10 and 7, respectively. 

For all 1 .3 < z < 2.4 objects, we estimated the AB absolute 
magnitude M2000 using a synthetic 1900 A<A<2100A box 
filter derived from the observed V and / magnitudes (Fig. 1). 
To estimate M2000 for the galaxies in our sample we use the 
following empirical formula: 



where: 



M 2000 = y-2/(1.72)±/(z)±5.058-2.91z 
-0.114(V-/)(z-1.72) 



/(z) =25±51ogD L (z)-2.51og(l + z) 



(1) 



(2) 



and Di(z) is the luminosity distance in Mpc of galaxies at 
redshift z- Eqs. 1 and 2 were derived by us considering that 
to first order approximation the observed V band matches the 
rest M2000 band at z = 1.72. The extra terms represent the 
mean ^-correction away from z = 1 .72 and a small correc- 
tion for its dependence on spectral energy distribution (SED) 
shape (parameterized by V -I observed). The coefficients 
were derived by empirically fitting the observed colors to 
1.3 <z < 2 template spectra (Chen et al. 2002) ranging from 
E/S0 to Starburst SEDs; the formula was found to predict 
■^2000 to < 0.2 mag accuracy over this entire redshift range for 
all SEDs. We emphasize that it is robust, model independent 
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(i.e. insensitive to SED shape) in this redshift range simply 
because the optical photometry samples the rest-frame UV. 

Among the 34 galaxies with secure redshifts, we selected 
the brightest subsample with M2000 < -19.65 (Fig. 1). This 
subsample contains 17 galaxies, but three of them were not 
considered because the spectra are contaminated by second 
order lines, and a fourth because it has partial spectral cover- 
age in the 2200-2900 A range. The remaining 13 galaxies 
(Table 1) cover the range 1.34 < z < 1.97 (mean and median 
redshift z = 1.60 and z = 1.56). In Table 1, we report redshifts 
independently derived from the ISM absorption lines and the 
[O II] emission line (when observed). The mean velocity dif- 
ference between the [O II] emission and ISM absorption is 
Av = 93 with 197 km s -1 dispersion. The observed I-K col- 
ors (in the range 2.2 < I-K < 3.5) are consistent with predic- 
tions derived from a suite of spectral synthesis models for Irr- 
or Scd-type galaxies. This sample represents 25% or 38% 
of the total number of GDDS galaxies in 1.3 < z < 2.4, de- 
pending whether the objects with tentative spectroscopic or 
photometric-only redshifts are included in the counting or 
not. 

4. THE [O II] EMISSION AND STAR FORMATION RATES 

The [O II] A3727 emission line is detected in 11 of the 13 
galaxies selected (Fig. 2). For the two highest redshift galax- 
ies the [O II] wavelength is out of the observed spectral range. 
The presence of the [O II] nebular emission is another indica- 
tion of star formation activity in these galaxies, as it is asso- 
ciated with H II regions. No other emission line is detected. 
The [O II] flux has been measured in the spectrum of each 
galaxy and reported in Table 1. As fluxes in the red end of 
the spectra are roughly consistent with the aperture corrected 
photometry, no aperture correction has been applied. Errors 
are derived using photon counting statistics only. Our spectral 
sensitivity is sufficiently good to detect very low fluxes, down 
to - 2 x 10" 17 erg s" 1 cm" 2 (~ 3cr) at A > 8600 A. Observed 
fluxes, in the range (3-26) x 10~ 17 erg s -1 cm -2 and absolute 
luminosities (5 - 42) x 10 41 erg s -1 are reported in Table 1. 

The star formation rate (SFR) in local H II regions is 
generally directly derived from optical Balmer lines (Gal- 
lagher, Hunter, & Bushouse 1989; Kennicutt 1992; Gallego 
et al. 1995). However, at high redshift, Balmer lines are red- 
shifted to the IR, making the SFR measurement harder. If 
Balmer lines are not observed, the [O II] emission line can be 
used instead, since it is detectable in the optical up to z ~ 1 .6. 
The problem is that [O II] gives a much less direct estimate 
of the SFRs, because it is more affected by dust than Balmer 
lines (for which dust correction can be easily determined), and 
because it depends on the electron temperature and metallic- 
ity of the H II region. In local galaxies, Jansen, Franx, & 
Fabricant (2001) found a correlation between [O II]/Ha ratio 
and the galaxy absolute magnitude that allows to estimate the 
SFR if only the [O II] emission line and the Mb are known. 
Their Fig. 4a shows that for bright objects (Mb < -20.5, as 
our galaxies are) Ha is about 5- 10 times brighter than [O II]. 
The SFR obtained with this method is uncertain by a factor of 
5. 

However, this relationship might not be valid at high red- 
shift and the SFR may be overestimated by up to a factor of 
10 (Jansen et al. 2001; Hicks et al. 2002). Glazebrook et al. 
(1999) have compared Ha and [O II] luminosities in a sample 
of 13 z ~ 1 field galaxies. The median value in the sample 
is Lhq ~2x L[ 011 j ; in this relation the dust correction is not 



applied. 

If we use the prescription (Kennicutt 1998): 



SFR HQ (M yr- 1 )~ 



1.26 x 10 41 ergs" 



(3) 



that gives the SFR for a given Ha luminosity, and assume a 
dust correction (see §8.2) for the Ha flux of Au a = 0.5 mag- 
nitudes (also adopted by Glazebrook et al. 1999), we obtain: 



SFR IniI1 (M Q yr _1 ) ~ r 

[on] 1 oy ) o.39 x 10 41 erg s" 1 



(4) 



where the [O II] luminosity is the one observed, i.e. not cor- 
rected for dust. From this, we derive SFRs in out sample in 
the range 13- 106 M yr -1 (Table 1), with mean and median 
value of <~ 40 M yr -1 . These values are compared in §8.3, 
with the SFRs derived using UV luminosities. 

High SFR values (derived from Balmer lines) are not new to 
high redshift star-forming galaxies. For instance, Glazebrook 
et al. (1999) derived 20-60 M yr" 1 in their z ~ 1 sample. 
From UV fluxes, Daddi et al. (2003) derived much higher 
SFRs in the K20 1.7 < z < 2.3 galaxies, with SFRs ranging 
100-500 M yr" 1 . High SFRs are also found in 2.3 < z < 3.3 
LBGs, for which infrared spectra indicate SFR ~ 20-200 M 
yr" 1 (Pettini et al. 1998; Kobulnicky & Koo 2000). 

5. CO-ADDITION 

In order to improve the signal-to-noise ratio in the 2200 - 
2900 A spectral range, and to enable studies of the global 
properties of the galaxy sample, the 13 selected GDDS spectra 
have been shifted to the rest frame and combined. Both lin- 
early and optimally extracted spectra were initially used and 
two separate composites have been obtained. To combine the 
data in a suitable way and maximize the signal-to-noise ratio, 
we have used the following procedure. The spectra were first 
normalized to one; this allowed us to analyze in the compos- 
ite spectrum the EWs of lines of the same ion in a consistent 
way. The continuum shape was determined by interpolating 
regions free of spectral features. Then we created for each 
spectrum a mask spectrum containing "0" or "1", according 
to whether the region or pixel can be used or not. Generally, 
the spectra are relatively free of masked pixels, however small 
regions in some spectra are contaminated by strong second or- 
der refraction. Occasionally other pixels were masked due to 
strong sky emission line residuals or CCD defects. The total 
number of masked pixels in the A = 2200-3000 A spectral 
range is 3.8% of the total. 

The observed spectra were shifted to the rest-frame (using 
ISM redshifts) and linearly averaged to form a mean spectrum 
(i.e. each galaxy was weighted the same in the composite). 
The composite noise spectrum is derived consistently. We fi- 
nally re-normalized the composite spectrum to adjust for mis- 
placements of the continuum in the individual spectra. The 
uncertainty in the continuum of the composite is estimated to 
be <~ 3-5%, from the higher to the lower signal regions, and 
is consistent with the level of the noise per resolution element. 

The composite optimal spectrum is shown in Fig. 3 (lower 
spectrum). For comparison, we also show (shifted upwards 
along the y-axis) the normalized composite of 14 spectra of 
local starburst galaxies (Tremonti et al. 2003) obtained with 
the low resolution Faint Object Spectrograph (FOS) on board 
of the Hubble Space Telescope (HST). The GDDS mean 
FWHM in the rest frame is 6.1 A, and is estimated by consid- 
ering the observed FWHM divided by (1 + (z)), where (z) — 
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FIG. 2. — [O II] emission for 11 GDDS galaxies. Smooth lines are the noise spectra per pixel (~ 3.5 A in the observed frame). Also reported is the central 
observed [O II] emission line wavelength A c . Note that the continuum is not subtracted from the spectra, but subtracted when measuring line fluxes. The flux at 
A ~ 9960 A in object GDDS-02-2264 goes below zero due to red-end contamination. 



1 .60 is the mean redshift of the sample. For the composite we 
have chosen a small pixel size of 0.7 A (Av ~ 80 km s" 1 at 
2600 A), i.e. ~ 2 times smaller than the real mean pixel size 
(1.39 A). This allows us to better identify small features in 
the spectrum and more accurately measure line EWs. The 
resolution element is extended over ~ 9 pixels. The compos- 
ite optimal noise spectrum per 0.7 A pixel is also shown as a 
dotted line in Fig. 3. 

As a consistency check, the EWs measured from the opti- 
mal and linear composite spectra were compared. The agree- 



ment is very good and no correlation with, e.g., wavelength is 
present. Since the signal-to-noise of the optimally extracted 
composite is 5-15% higher than the linearly extracted one, 
the optimal composite is used for all subsequent analysis. 

The mean S/N per 0.7 A pixel in regions free of strong fea- 
tures has been estimated using the noise spectrum. This is S/N 
= 9.1 at A = 2210-2330 A, S/N = 12.3 at A = 2395-2570 
A, S/N = 14.2 at A = 2615-2785 A, S/N = 11.8 at A = 
2815-2990 A. The S/N measured with the noise spectrum 
rescales with the square-root of the pixel size chosen. As a 
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FIG. 3. — Composite spectrum of 13 GDDS galaxies with ISM absorption lines (lower spectrum). The noise per 0.7 A pixel (dotted spectrum) is also shown. 
The sample covers the redshift range 1.34 < z < 1.97 (mean and median redshift z =1.60 and 1.56, respectively). Absorption features are also marked. As 
a reference, we show the composite spectrum of 14 local starburst dwarf galaxies observed with HST/FOS (Tremonti et al. 2003), shifted along the y-axis for 
comparison. 



sanity check, we also measured the standard deviation directly 
in the signal spectrum. We obtain S/N = 1 1.2, 13.6, 18.3, and 
13.9, in the same regions as before, i.e. about 1.2 times higher 
than that obtained from the noise spectrum. This is close to 
the expected value, y/2 ~ 1 .4, given that the chosen pixel size 
(0.7 A) is ~ 2 smaller than the real pixel size (1.39 A). 

6. HEAVY ELEMENT ENRICHMENT 

The ultimate goal of the present analysis is to estimate the 
heavy element enrichment of galaxies in the GDDS compos- 
ite spectrum. In the observed range, we have clearly detected 
Fe II, Mg II and Mg I absorption lines (Fig. 3). Besides these, 
we have also included in the analysis the very weak Mn II 
triplet around 2600 A. All EWs with errors in the GDDS com- 
posite are reported in the third column of Table 2. The 4a EW 
upper limits for two Fe II non-detections are also estimated 
using the noise spectrum, by assuming that lines are extended 
over 6 A (or Av ~ 800 km s" 1 at A - 2250 A). Values in 
brackets are EWs obta ined by extrapolating from estimated 
column densities, see § 16.11 To check the line de-blending in 
the Fe II AA2374,2382 and Mg II AA2796,2803 doublets, we 
compared the measured EWs in the spectrum with the EWs 
given by the fit of two blended Gaussians. For the Fe II dou- 
blet we measure a ratio between the two lines of 0.59, to be 
compared with 0.64 given by the ratio of the two Gaussians. 
For the Mg II doublet we obtain 1.05, to be compared with 
1.03 given by the two Gaussians. This shows that the con- 
tamination is properly accounted for. The contamination of 
the Fe II A2600 line by Mn II is treated separately (see next 
section). 

The UV continuum of star-forming galaxies (like our se- 
lected sample) is dominated by O and B stars. In this case, the 



stellar contamination to the ISM lines in the 2200 - 3000 A in- 
terval is negligible (Tremonti et al. 2003). Moreover, when 
stellar wind is not important (this is the case for lines in the 
2000-3000 A region) the strength of the photospheric lines is 
significantly weaker when compared to ISM lines (de Mello, 
Leitherer, & Heckman 2000). 

The lines detected in the GDDS composite are strong, but 
narrow in velocity space (Av < 800 km s -1 ). Assuming that 
these interstellar lines originate in gas viewed against the con- 
tinuum provided by the integrated light of O and B stars in the 
galaxies, we can measure a sort of "average" column densi- 
ties of the ions. Two approaches have been followed. In the 
first case we used EWs as given in Table 2 with the curve of 
growth analysis (COG, Spitzer 1978). In the second, we si- 
multaneously fit absorption profiles to all lines in the spectrum 
arising from transitions of the same ion. The advantage of fit- 
ting line profiles is that line blending is automatically taken 
into account. 

6.1. Column densities 

The rest-frame EWs of the ISM lines detected in the GDDS 
composite galaxy are very high, indicating that the gas is 
heavily enriched by metals. As the majority of absorption 
lines are very strong (W r > 1 A), they deviate from the lin- 
ear part of the COG and it is not possible to use the linear 
approximation to determine column densities. Moreover, the 
low resolution (FWHM ~ 700 km s" 1 at A - 2600 A) does 
not allow application of the apparent optical depth method 
(Savage & Sembach 1991). However, even with the relatively 
large uncertainties, the general COG analysis can still give 
significant results since so many lines are present. Given that 
several galaxy spectra are combined together, it is not pos- 
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Table 2 

Rest-frame equivalent widths 



Line 


fx 




W r (A) 




GDDS 


cB58 a 


LBGs 


Si II A 1260 


1.0070 




(2.54 ±0.07) 


< 1.7 


Si II A1526 


0.1270 




2.59 ±0.03 


1.72±0.18 


Si II A 1808 


0.00218 




0.53 ±0.02 


(0.04-0.28) 


Fell A 1608 


0.058 


(1.62±0.12) 


1.30±0.04 


0.91 ±0.15 


Fe II A2249 


0.00182 


< 1.2 


(0.15±0.01) 


(0.048-0.13) 


Fe II A2260 


0.00244 


< 0.9 


(0.19±0.01) 


(0.065-0.18) 


Fe II A2344 


0.114 


2.90 ±0.27 


2.99 ±0.04 




Fe II A2374 


0.0313 


2.28 ±0.23 


1.94 c 


(0.9-1.4) 


Fe II A2382 


0.320 


3.57±0.31 


3.41 d 




Fe II A2586 


0.0691 


2.96±0.17 


(2.88 ±0.06) 




Fe II A2600 


0.239 


3.94 ±0.27 


(3.95 ±0.10) 




Mn II A2576 


0.3508 


0.56 ±0.1 3 


0.38±0.10 




Mn II A2594 


0.2710 


(0.45 ±0.16) 


(0.30 ±0.09) 




Mn II A2606 


0.1927 


(0.33 ±0.12) 


(0.22 ±0.06) 




Mg II A2796 


0.6123 


4.10±0.I9 






Mg II A2803 


0.3054 


3.98±0.16 






Mg I A2852 


1.8100 


0.91 ±0.19 







NOTE. - Oscillator strengths fx are taken from a recent compilation by 
Prochaska et al. (2001). EWs between brackets are extrapolated from mea- 
sured column densities (see text). EW upper limits for Fe II in GDDS com- 
posite are 4<r. 

"EWs taken from Pettini et al. (2002). 

b EWs taken from Shapley et al. (2003). 

c The lcr error reported by Pettini et al. (2002) is > 0.07 A (line is affected 
by bad pixels). 

d The lcr error reported by Pettini et al. (2002) is > 0.05 A (line is affected 
by bad pixels). 
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FIG. 4. — Best fit curve of growth for Fe II in the GDDS composite spec- 
trum. The solid and dashed lines are for the best fit Doppler parameter and 
±lo\ respectively. The other ion column densities are calculated assuming 
the best fit b value found for Fe II. The straight line shows the linear approx- 
imation to the curve of growth (b = +°°). 
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FIG. 5. — Line absorption fitting profile of Fe II for GDDS composite. 
The smooth line is the best fit with logAfceii = 15.52+g;}5 and b = 1 15 ± 10 
km s . Shown at the bottom are the best fit residuals. 



sible to separate the different components of very complex 
absorption features and therefore they appear as single lines 
(Fig. 3). Hence, the widths of these features do not provide 
useful information on the kinematics or thermal conditions of 
the gas in the galaxies and the broadening is a sort of "effec- 
tive" Doppler parameter, being the result of the superposition 
of many absorbing clouds. The COG analysis with a single 
component approximation was used for QSO-DLA studies 
before high resolution spectroscopy of QSOs became possi- 
ble (see for instance Blades et al. 1982; Turnshek et al. 1989; 
Steidel 1990). In fact, Jenkins (1986) showed that the COG 
technique applied to complex features nearly always gives a 
reasonable answer (the simulated-to-true column density ra- 
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tio rarely goes below 0.8) even if different lines have very 
different saturation levels or Doppler parameters. The general 
results have been confirmed by Savage & Sembach (1991), 
who in addition, show that the COG works best if the gas ab- 
sorption is not "bimodal" - e.g., that very large absorptions 
are not combined with narrow ones. Since our composite 
spectrum is derived from many galaxies, each with many ab- 
sorption clouds, a bimodal distribution of absorption lines is 
very unlikely. 

In Fig. 4 we show the COG for the absorption lines of the 
GDDS composite spectrum. The fit to the data was deter- 
mined by x 2 minimization of the COG. The resulting column 
densities are reported in Table 3. To determine the Fe II col- 
umn density, we have used 5 strong detections, and iterated 
the results with those for Mn II, because two of the three weak 
Mn II lines are blended with the Fe II A2600 line. The weakest 
Fe II line detected, Fe II A2374, gives logA^eii > 15.2, using 
the linear approximation for the COG. The non-detection of 
Fe II A2260 line sets an upper limit of logA^i! < 16. These 
two limits already give a good, robust estimate of the Fe II col- 
umn density. The 5 Fe II EWs give a best fit column density 
of logAfpeii = 15.54^q j3, and an effective Doppler parame- 
ter b = 111 =b 11 km s . From this, we extrapolate the EW 
of the Fe II A1608 line to be W r = 1.62 ±0.12 A. The re- 
duced x 2 value obtained from the best fit is 0.34, indicating 
that the errors are probably not underestimated. The Mn II 
A2576 line is isolated and marginally detected with 4.4er sig- 
nificance (Table 2). Moreover, the presence of the Mn II ab- 
sorption is supported by the large EW of the Fe II A2600 line, 
indicating the presence of Mn II AA2594,2606 there too. The 
column density derived from the Mn II A2576 line is in the 
range 13.25 < logA^vmii < 13.67, assuming b > 90 km s , 
with very little saturation. Taking the same b value as found 
for Fe II, we find a best fit logjV Mn ii = 13.49 ± 0.13. From 
this, we derive W r = 0.45 ±0.12 A and 0.33 ±0.10 A, for 
Mn II A2594 and Mn II A2606, respectively. These two val- 
ues have been subtracted from the Fe II A2600 EW and the 
error propagated. The weak Mg I A2852 line gives logA^Mgi = 
12.93 ± 0.13 for b = 111 km s" 1 . The Mg II doublet is 
heavily saturated, namely logA^Mgii > 14.3. We estimated 
log/Yjy[gii = 15-26+q 24 for a Doppler parameter b — 1 1 1 ± 1 1 
km s -1 (as for the Fe II best fit). 

As a further test on the Fe II column density measurement, 
we have performed simultaneous line profile fitting of the 
5 lines. The result with one component gives logA/peii = 
15.52^ }^ and b = 1 15 ± 10 km s _1 (Fig. 5), very similar to 
that found with the COG method (log7Y FeII = 15.54+|] ^), with 
a reduced x 2 — 0.55. This shows consistency between the line 
profiles and the column density derived using the COG. The 
velocity shift between the fit and the rest frame is consistent 
with zero. 

We finally considered another source of uncertainty: the 
ISM covering factor. So far we assumed that the covering 
factor is one, that is the gas is uniformly distributed in front 
of the emitting sources, and photons escaping the absorbing 
gas or scattered are not important. However, in case of a non- 
unity covering factor, the residual flux not absorbed by the 
ISM would cause a systematic underestimate of the column 
densities, if not properly taken into account. From the depth 
of the Mgll absorption, we have estimated that the covering 
factor is at least 70%. We recalculated the Fell column den- 
sity, assuming this very conservative lower limit, and obtained 



Afceii ~ 60% higher than in the case of unity covering factor. 
In log scale the difference is only 0.2 dex: log/A^eii = 15.75 
instead of logAfp e n = 15.54; this is within the upper error 
estimated in our first measurement (0.23 dex). Shapley et 
al. (2003) discuss the non-unity of the covering factor in 
LBGs. However, Giallongo et al. (2002) found that in two 
LBGs the escaping fraction of ionizing photons (A < 912 A) 
is < 16%, whereas Malkan, Webb & Konopacky (2003) found 
that the same quantity in the UV images of 1 1 bright blue 
1.1 < z < 1.4 galaxies is less than 6%. Both results evoke 
some caution on this issue. 

7. METAL COLUMN DENSITIES IN HIGH REDSHIFT ISM 

When the gas is mostly neutral (A/hi ~ Nh), the column den- 
sity of an element can be approximated with the column den- 
sity of the ion with ionization potential above the hydrogen 
ionization potential (13.6 eV). For instance, iron, magnesium 
and manganese are mainly in the Fe + , Mg + , and Mn + state. In 
the GDDS composite, the very high EWs of Fe II and Mg II 
indicates that the gas is nearly neutral and the Np s u ~ A/p e , 
A^Mgii ~ A^Mg, or A^Mnii ~ A^Mn approximation can be applied. 
The low ionization of the gas is confirmed by the presence of a 
relatively strong Mg I absorption (logA^Mgi ~ 12.9), but with 
column density much lower than the Mg II column density 
aogJV Mg n~ 15-15.6). 

Considering column densities of these ions in the neu- 
tral ISM instead of metallicities (element-to-hydrogen relative 
abundance) has the advantage that we can directly compare 
results with other classes of objects at various redshifts with 
similar measurements, without knowing the H I content. In- 
deed, H I absorption in GDDS spectra cannot be detected for 
lack of spectral coverage. In Table 3 we compare column den- 
sities of Fe II, Mn II, Si II, observed in the GDDS composite, 
and other high redshift and local galaxies. 

Heavy element column densities in the neutral ISM of high 
redshift galaxies have been measured in great detail for more 
than 100 QSO-DLAs in the interval 0.0 < z < 4.5 (see for in- 
stance Lu et al. 1996; Pettini et al. 1997; Prochaska & Wolfe 
1999; Prochaska et al. 2003a). For our purpose, we have used 
our own large compilation of DLAs 18 . The mean column den- 
sities for Fe II, Mn II, and Si II with dispersions are given in 
Table 3. We have divided the sample in subsamples for each 
ion. The number of DLAs considered in each subsample is 
also given in the first column (between brackets). The col- 
umn density distributions for Fe II and Mn II as a function of 
redshift are shown in Figs. 7 and 8. 

The heavy element column densities in the ISM of GRB- 
DLAs have been studied in only 4 objects (Savaglio et 
al. 2003, Vreeswijk et al. 2003). In a fifth event, GRB 020813 
(z = 1 .255), equivalent widths of many UV lines have been re- 
ported by Barth et al. (2003), from which we have estimated 
column densities using the COG (Table 3). 

We also include in our analysis, the ISM lines measured in 
star-forming galaxies in the local Universe observed with the 
Far Ultraviolet Spectroscopic Explorer (FUSE). The sensitiv- 
ity and resolution (FWHM ~ 30 km s" 1 ) of FUSE is good 
enough for a detailed analysis of nearby starburst galaxies. 
In our study we considered metal column densities obtained 
for NGC 1705 (Heckman et al. 2001) and I Zw 18 (Aloisi et 
al. 2003). For comparison, we have also included measure- 
ments in the diffuse and dense ISM of the Milky Way and 

18 For the largest public compilation of DLA column densities, visit 
J. Prochaska's homepage at http://www.ucolick.org/~xavier/ESI 
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FIG. 6. — Best fit curve of growth for Fe II in the cB58 ISM. The solid and 
dashed lines are for the best fit Doppler parameter and ±1ct, respectively. 
The Mn II column density is calculated assuming the best fit b value found 
for Fe II. The best fit for Si II is calculated independently and is not shown 
for clarity, but this is very similar to the Fe II best fit. The straight line shows 
the linear approximation to the curve of growth (b = +°°). 



the Magellanic Clouds (Sembach & Savage 1996; Spitzer & 
Fitzpatrick 1995; Welty et al. 1997, 1999). 

The neutral ISM metal enrichment of LBGs has been dis- 
cussed by Shapley et al. (2003) for a composite of 8 1 1 spectra, 
and by Pettini et al. (2002) for the bright LBG galaxy cB58 
at z = 2.72. In the following sections, we use their EWs and 
estimate heavy element column densities using the COG. 

7.1. The spectrum of the LBG MS 1512-cB58 

Column densities in the LBG MS 1512-cB58 (z = 2.72) 
have been determined by Pettini et al. (2002) using the ap- 
parent optical depth (Savage & Sembach 1991) measured on 
the less saturated lines. For a comparison of the COG with 
this method, we have redone the analysis for cB58 using the 
COG. The EWs used are given in the fourth column of Ta- 
ble 2 and are the ones reported by Pettini et al. (2002). The 
COG for cB58 is shown in Fig. 6 and the column densities are 
listed in Table 3. The EWs of two of the four Fe II lines, Fe II 
AA2374,2382, have errors larger than 0.07 and 0.05 A for 
contamination by bad pixels (Pettini et al. 2002), for which 
we arbitrarily use a = 0.15 A. 

The advantage of using the COG is that all detected Fe II 
lines can be used, regardless of the saturation level. The ap- 
parent optical depth method, as Pettini et al. (2002) points 
out, can still underestimate column densities due to satura- 
tion. This effect is also noted in the ISM of local dwarf 
galaxy I Zw 18 by Aloisi et al. (2003), who made a com- 
parison between the line fitting of Voigt profiles and the ap- 
parent optical depth method. The COG analysis for cB58 
is shown in Fig. 6. The best fit to the four Fe II lines gives 
logA? F eii 15.25±g$ and b = 130 ±4 km s" 1 ( X 2 = 3.9), con- 
sistent with logJVpeii — 15.20 found by Pettini et al. (2002). 



Errors on column densities derived by Pettini et al. (2002) are 
reported to be generally less than 0.2 dex. 

The Mn II A2576 line is very close to the linear part of 
the COG, therefore the column density is relatively well con- 
strained in the range log_/V M nii = 13.12 - 13.44 for b > 90 
km s . If we adopt b = 130 km s -1 , as given by Fe II, we 
obtain logA^vinii = 13.29 ±0.13, consistent with logA^vinii — 
13.33 found by Pettini et al. (2002). From this column density 
and Doppler parameter, we derived the expected EWs of the 
Mn II AA2594,2607 doublet (Table 2). 

For completeness, we also estimate the Si II column den- 
sity. Si II A 1260 is weakly contaminated by S II A 1259 and 
Fe II A 1260 (W r ~ 1 A), therefore we do not use it. The other 
two detected lines have very different oscillator strengths. 
The Si II A 1808 line is weak and little saturated, and gives 
a column density in the range 15.92 < logA^sm < 16.05 for 
b > 70 km s _1 . From Si II AA1526, 1808 together, we de- 
rive logA^sm = 15.99+0 02 and ^ = 136 ± 2 km s" 1 , consistent 
with logA^sin ~ 15.99 given by Pettini et al. (2002). We also 
note that the effective Doppler parameter for Fe II and Si II 
have been found independently to be very similar (bp e n = 130 
km s -1 and bsm = 136 km s -1 ). 

We found an excellent agreement between the results ob- 
tained by Pettini et al. (2002) using the apparent optical depth 
method and ours obtained with the COG (Table 3), supporting 
the robustness of our results for GDDS galaxies. 

7.2. The LBG composite spectrum 

Shapley et al. (2003) have used 811 LBG spectra in the red- 
shift range 2 < z < 3.8 to construct a UV (A = 950- 1850 
A) composite spectrum. Their composite is rather complex, 
showing stellar features, low and high ionization absorption, 
and nebular emission lines. EWs of ISM lines are reported, 
but column densities are not derived because lines are heavily 
saturated. They also subdivided the sample into 4 subsam- 
ples according to the Lya equivalent width, in order to find 
possible correlation with observed parameters. Here we only 
consider the EWs of the total composite. Although the few 
ISM lines are saturated (W r > 0.9 A), inspecting them with 
the COG can provide interesting, even if qualitative, results on 
the heavy element enrichment in the neutral ISM. We assume 
that the covering factor of the ISM in front of the UV sources 
is one. If the covering factor is less than one, column den- 
sities derived using the COG would be systematically lower 
than real values. 

Only Fe II A 1608 is measured for Fe II, but because this line 
is moderately saturated, it can be used to constraint the col- 
umn density. We derive 14.76 < logA% e n < 15.21, assuming 
an effective Doppler parameter b > 90 km s" 1 (a reasonable 
assumption given that 811 galaxy spectra are coadded). From 
this, we derive the expected EWs for other Fe II transitions 
(between brackets in Table 2). The Si II A1304 line is con- 
taminated by O I A 1302, while Si II A 1260 is contaminated by 
S II A 1259 and Fe II A 1260. To determine the Si II column 
density, we used these two EWs as upper limits and in addi- 
tion the Si II A 1526 measurement. Assuming ^ > 90 km s , 
we get logA^sm = 14.8- 15.7, from which we expect W,-(Si II 
A1808) = 0.04-0.28 A (Table 2). 

A comparison between the GDDS composite and the LBG 
composite can be made by considering the Fe II A 1608 line. 
Given the good accuracy of our Fe II column density measure- 
ment, we can estimate the expected EW for the Fe II A 1608 
line (that is not covered by our data), and get W r — 1 .6 ± 0.2 
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Table 3 

Interstellar column densities 



logN (cm 2 ) 



Galaxy 


redshift 


Fell 


Mn II 


Si II 


HI 


Ref 


V jULtj y id) 


i £ft 
1 .ou 


D ' M -0.13 


1 3 AQ -4- ft 1 3 
13. zt U. 1 j 






i 
i 


MW-MC 


0.00 


(14.89) 


(13.00) 


(15.53) 




2,3,4 


IZw 18 


0.0026 


15.09 ±0.06 




14.81 ±0.07 


21.34±0.11 


5 


NGC 1705 


0.0021 


14.54 ±0.10 




14.61 ±0.14 


20.2±0.2 


6 


GRB 990123 


1.6004 


14.78^12 








7 


GRB 000926" 


2.038 






16.47+jU" 


-21.3 


7 


GRB 010222 


1.475 




13 6l +aos 


16.09 ±0.05 




7 


GRB 0208 13 b 


1.255 


. 5 59+0.02 


13.63 ±0.02 


16.28 ±0.03 




1 


GRB 030323 


3.3714 


15.93 ±0.08 






21.90±0.07 


8 


cB58 c 


2.72 


, 5 25+ ' 04 


13.29±0.13 


1 5 on+0.01 
13 - yy -0.02 


20.85 ±0.09 


9 


LBGs (811) c 


~3 


14.8-15.2 




14.8-15.7 




10 


QSO-DLAs (120) d 


0.56-4.47 


14.39 (0.62) e 






20.60 (0.51) e 


11 


QSO-DLAs (21) d 


0.56-2.78 




12.68 (0.44) e 




20.68 (0.56) e 


11 


QSO-DLAs (105) d 


0.68-4.47 






14.71 (0.70) c 


20.52 (0.50) e 


11 



References. - (1) This work; (2) Sembach & Savage (1996); (3) Spitzer & Fitzpatrick (1995); 
(4) Welty et al. (1997, 1999); (5) Aloisi et al. (2003); (6) Heckman et al. (2001); (7) Savaglio et al. 
(2003); (8) Vreeswijk et al. 2003; (9) Pettini et al. (2002); (10) Shapley et al. (2003); (11) Obtained 
from a compilation of QSO-DLAs 

a H I column density estimated by Fynbo et al. 2002. 

b Line EWs are taken from Barth et al. 2003. 

°We have estimated column densities using the COG. 

d In brackets is the number of QSO-DLAs used. 

e In brackets is the dispersion in the sample. 



A (Table 2). This is 1 .8 times higher than the EW of the same 
line measured by Shapley et al. (2003) in the LBG composite 
(W r = 0.91 ±0.15 A), but it is similar to that of the LBG sub- 
set (1/3 of the galaxies) that does not display Lya emission 
(W r = 1.57 ±0.21 A). 

8. DISCUSSION 

Various methods have been traditionally used to derive 
metallicities in high redshift galaxies using both absorption 
and emission lines. The Lick indices method (Faber et al. 
1985), recently applied on SDSS 0.15 < z < 0.50 galaxies 
(Eisenstein et al. 2003), compares absorption line strengths 
in stars with models (Vazdekis 1999) and is mostly sensitive 
to old stellar population systems. As Eisenstein et al. (2003) 
point out, the models are limited in the range of metallicities 
and can suffer from the age-metallicity degeneracy. 

Nebular emission in the ionized ISM is, on the other hand, 
more appropriate when dealing with star-forming galax- 
ies. The metallicity, measured using [O II] A3727, [O III] 
AA4959,5007, and Balmer emission lines of 56 galaxies at 
0.26 < z < 0.86, ranges from half to two times solar (Kob- 
ulnicky et al. 2003). In a similar study (66 galaxies at 
0.5 < z < 0.9) Lilly, Carollo & Stockman (2003) find that 
most galaxies in their sample have Z <~ Z©. This method re- 
lies on a good calibration of metallicity vs. line strengths de- 
rived from samples in the local Universe, and requires also a 
reliable estimate of the stellar Balmer absorption line contam- 
ination. As Kobulnicky et al. (2003) mention, it is useful in 
a statistical sense when using large data sets. Moreover, this 
method requires IR spectroscopy for galaxies at z > 0.8. 



In our first analysis of GDDS data, we used UV absorption 
lines associated with the neutral ISM to evaluate heavy ele- 
ment enrichment in a subset of galaxies in 1.34 < z < 1.97. 
This sample represents a joint set of the most massive galax- 
ies at z > 1.3 (from the deep K selection) together with those 
with the strongest rest-frame UV emission (from the M2000 
selection). Thus we are selecting massive starburst galax- 
ies; preliminary stellar mass estimates (Glazebrook et al. in 
preparation) indicate that they have masses M ^ 10 10 M Q and 
their descendants must correspond to galaxies of at least this 
mass today. This selection is quite similar to that of Lilly 
et al. (2003) at z ~ 0.7 whose primary /-band selection pulls 
out massive galaxies and whose emission line cut pulls out a 
high-SFR subset. 

8.1. Metal content at high redshift 

One of the advantages of using ISM absorption lines associ- 
ated with singly ionized elements in galaxies is that no model- 
dependent assumptions or calibrations, dust extinction or ion- 
ization corrections are important when deriving the heavy el- 
ement enrichment. 

In Fig. 7 we show the Fe II column density measured in the 
ISM of different galaxies, as a function of redshift. In our 
analysis we unambiguously find that the metal enrichment of 
the GDDS galaxies is high, with log/V Fe ii = 15.54^. Eigh- 
teen QSO-DLAs have measured Fe II column density in the 
same redshift range (1.34 <z< 1.97). In this sample the high- 
est column density of Fe II is logA/peii — 15.36, i.e. ~ 1.5 
times lower than in the GDDS composite. If we extend the 
range to 1 .0 < z < 2.0, we find 24 QSO-DLAs, and 2 systems 
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FIG. 7. — Fe II column densities as a function of redshifts in differ- 
ent galaxies. QSO-DLAs are divided in three subsamples where small, 
medium and big squares indicate logWHi < 20.45, 20.45 < logWui < 20.80, 
and log^Hi > 20.80 systems, respectively. Errors for logAfcen in QSO- 
DLAs is always < 0.2 dex, and on average is 0.05 dex. Errors for lo- 
cal ISM column densities, when available, are typically less than 0.1 dex. 
For logA'peii in the LBGs, we have taken the middle point in the interval 
logAfa-u = 14.76— 15.21. Errors on the other column densities are given in 
TableTy 



14 



□ 



13 - 



o 



12 



• GDDS 

xMW * MC ISM 
□ □ □ QSO-DLAs 
» GRB-DLAs 

CMS 1512-cB58 



3 



2 

redshift 



FIG. 8. — As in Fig. 7, but for Mn II column densities. The error for 
logWiymil in QSO-DLAs is on average 0.07 dex. Errors for local ISM column 
densities, when available, are typically less than 0. 1 dex. 



with comparable column density (8% of the total, or 15% if 
only logA^Hi > 20.3 systems are considered). This is still a 
significant difference with the GDDS composite. On the other 
hand, high values of Fe II column density are found in other 
star-forming galaxies: the GRB host galaxies, and the LBGs 
cB58atz = 2.72(Table3). 

We note that when the background source is external to the 
galaxy (as for QSO-DLAs) the column density of an ion is 
higher than when the background source is inside the galaxy 
(as for the GDDS composite). In the former case, the line 
of sight crosses the entire galaxy, whereas in the latter only 
the gas in front of the stars is probed. The simplest situation 
would be that the stars are at the center, therefore the observed 
column is a factor of 2 less than what a background QSO 
would see. For some distribution of stars and gas, the cor- 
rection could be larger or smaller, but never downward. This 
would reinforce our findings of high Fe II column density in 
GDDS galaxies with respect to typical Fe II column densities 
in QSO-DLAs. 

Similarly, in Fig. 8 we show the Mn II column densities 
in the ISM of different galaxies. Even if the number of 
detections is much lower, the trend shown by Fe II is con- 
firmed. In the whole QSO-DLA sample, only two (out of 21, 
9.5%, or 13% if logA^i > 20.3 systems only are considered) 
have Mn II column density comparable to the GDDS compos- 
ite (logAWn = 13-49 ±0.13). At 1 <z< 2, all 9 QSO-DLAs 
with detected Mn II have logA^tnH < 13.1. Comparable val- 
ues are found, once again, in GRB host galaxies and cB58 
(Table 3). 

The Mg II doublet is very strong in the GDDS compos- 
ite (Fig. 3), but is highly saturated and the column density 
very uncertain (logA^Mgii = 15. 26+024). From this, we de- 
rive a magnesium-to-iron relative abundance 19 in the range 
-0.84 < [Mg/Fe] < 0.13. In the Galactic ISM, Fe is more 
depleted in dust than Mg, and [Mg/Fe] in the gas form is 
in the range 0.1 < [Mg/Fe] < 0.8. A smaller content of Fe 
with respect to Mg is also expected in the case of a-element 
enhancement (Mg is an a element) in metal poor systems. 
The low value of [Fe/Mg] found in the GDDS composite may 
also indicate that there is a residual saturation in the measured 
Mg II column density that is not taken into account. 

The Mg II absorption has been extensively studied along 
QSO sight lines (Rao & Turnshek 2000; Churchill, Vogt, & 
Charlton 2003; Ding et al. 2003). These absorbing clouds 
have been claimed to be associated with galaxy environment 
with logA^i k, 17.3. When the EW of the Mg II doublet is 
less than ~ 1 A, the absorbing gas originates in halo clouds, 
while for larger EWs, the system is very likely a DLAs (Rao 
& Turnshek 2000). In GRB afterglows, the Mg II absorption 
is, when detected, very strong and often used to identify the 
redshift, since it can be detected in the optical up to z ~ 2.2. 
In Fig. 9 we show the Mg II doublet EWs seen in QSO sight— 
lines, GRB afterglows, and GDDS composite, as a function of 
redshift. The difference between QSO sight-line studies and 
star-forming galaxies is apparent. 

Unfortunately for lack of wavelength coverage we have 
no information on the Lya absorption in our GDDS spec- 
tra. Therefore the H I column density and the metal content 
relative to hydrogen cannot be derived. A rough estimate of 
the H I column density can be obtained using the correlation 

19 We adopt the definition [X/Y] = logNx/Nr -log(X/Y) B , where 
(X/Y)q is the meteoritic solar abundance of element X with respect to el- 
ement Y as given by Grevesse & Sauval 1998. 
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between the SFR surface density and the H I surface den- 
sity, found by Kennicutt (1998) for local disk galaxies. Our 
slit aperture of 0.75 arcsec corresponds to a physical size (at 
z = 1 .6) of 6.4 kpc. The signal along the y direction has been 
extracted over 9.6 pixels (average value for the 13 galaxies) 
and this corresponds to 0.70 arcsec, or 5.9 kpc. If we assume 
that the detected flux comes from a 6.4 x 5.9 = 37.7 kpc 2 sur- 
face, and take the median SFR of the sample ~ 40 M Q yr" 1 
(derived from [O II] emission), we obtain a SFR surface den- 
sity: 



logS 



SFR : 



: 0.02 M yr" 1 kpc" 2 



(5) 



which corresponds to a H I surface density of log Shi ~ 1 .6 
M pc" 2 (±0.3 dex dispersion), or to a H I column den- 
sity \ogN m = 21 .7 ± 0.3. Changing the SFR by 20%, would 
change the H I of ~ 10%. Taking the Fe II column density de- 
termined with the COG (logA^eii = 15.54^j3), and assuming 
a moderate Fe dust depletion, Sp e = logA^F e .tot _ l°g-'VFe,gas = 
1 .0 dex (similar depletions are found in Galactic warm disk 
and halo+disk gas clouds; Savage & Sembach 1996), we ob- 
tain a metallicity Z/Z = 0.2+j] j. The error includes 20% and 
0.2 dex dispersion in the SFR and dust depletion, respectively, 
as well as the uncertainty in the H I and Fe II column densities. 

The median H I column density for 135 QSO-DLAs with 
logA^ m > 20.3 is logA^Hi ^ 20.7. About 1/4 of these have 
logA^ m > 21.0 and none has logJV H i > 21.85. For the LBG 
cB58 a logiVni — 20.9 have bee reported (Savaglio et al. 2002; 
Pettini et al 2002). To date, this is the only LBG for which 
the H I column density has been measured. High H I col- 
umn densities have also been estimated in GRB 000926 and 
GRB 030323 (Table 3). In addition, for GRB 000301C (z = 
2.040) and GRB 020124 (z = 3.198) a tentative logjV HI sa 21 .2 
and 21.7, respectively, have been reported (Jensen et al. 2001; 
Hjorth et al. 2003). If we assume for the GDDS composite a 
H I column density in the range logA^i ~ 21 .0-21 .7, and ap- 
ply the same dust correction as before (<5p e = 1 dex), we obtain 
Z/Z - 1-0.2. 

A high metallicity would not be terribly surprising, because 
even if QSO-DLAs dust-corrected metallicities are on aver- 
age ~ 1/10 solar in the redshift range 1 < z < 2 (Savaglio 
2000; Prochaska et al. 2003b), at redshift 0.3 < z < 0.9 higher 
metallicities (0.5 to 2 times solar) have been found in H II re- 
gions (Kobulnicki et al. 2003; Lilly et al. 2003). For cB58 at 
z = 2.72 Pettini et al. (2002) report Z/Z - 2/5. Heckman 
et al. (1998) have found that in the local Universe more lumi- 
nous starbursts have also higher metallicities and are dustier. 
Similarly, our subsample of GDDS galaxies is selected to be 
the most luminous in the sample, therefore they could also be 
the most metal rich. 

8.2. Dust depletion and extinction in the gas 

The relative abundances of elements with different deple- 
tion levels can be used to estimate the dust content in galax- 
ies. For example, Zn is typically hardly depleted in the ISM, 
while Fe is much more depleted. The zinc-to-iron relative 
abundance is on average [Zn/Fe] ~ 0.52, in 44 QSO-DLAs 
((z) = 2.0) with a \a dispersion of 0.28 dex. This is reminis- 
cent of the depletion pattern in the warm halo clouds of the 
MW, where [Zn/Fe] ~ 0.6 (Savage & Sembach 1996). 

In the GDDS composite, the Zn II doublet is off the wave- 
length range, whereas Fe II and Mn II were both detected. Un- 
fortunately, both Fe and Mn are strongly affected by the pres- 
ence of dust (Savage & Sembach 1996). We obtain a very un- 
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FIG. 9.— Equivalent widths of the Mgll AA2796.2803 doublet asafunc- 
tion of redshift in GDDS composite galaxy, QSO-Mg II absorbers, GRB 
afterglows and star-forming galaxy NGC 1705. QSO-Mg II absorbers at 
z < 1 .65 and z > 1 .65 are from Rao & Turnshek (2000) and Steidel & Sargent 
(1992), respectively. QSO-Mg II systems with total EW larger than W r = 1 
A (dashed line) are most likely QSO-DLAs. The four GRB detections are 
from, in order of increasing redshift, GRB 990712 (z = 0.433; Vreeswijk et 
al. 2001), GRB 020813 (z = 1.255; Barth et al. 2003), GRB 010222 (z = 
1.475; Jha et al. 2001; Masetti et al. 2001; Mirabal et al. 2002; Salamanca 
et al. 2001), GRB 990123 (z = 1.600; Kulkarni et al. 1999), GRB 990712 
(z = 1.619; Vreeswijk et al. 2001) and GRB 000926 (z = 2.038; Castro et 
al. 2003). The EW for NGC 1705 are from Vazquez et al. (in preparation). 
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FIG. 10. — Manganese-to-iron relative abundances as a function of Fe II 
column density in different high redshift galaxies and in the MW and MCs. 
The dashed lines indicate [Mn/Fe] found in, from top to bottom, cool disk, 
warm disk, warm disk+halo and warm halo clouds of the MW (Savage & 
Sembach 1996). 
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certain [Mn/Fe] = 0.07+[j ^ in the GDDS composite (Fig. 10). 
Future analysis with more data will allow to constrain bet- 
ter this measurement. In 28 QSO-DLAs with (z) = 1 .59, we 
found [Mn/Fe] = -0.02 with 0.22 dex dispersion. In the MW, 
the depletion patterns found by Savage & Sembach (1996) in- 
dicate [Mn/Fe] in the range 0.0-0.8, going from warm halo 
to cool disk clouds. A deficiency of Mn with respect to Fe 
has been found in low metallicity stars, due to "odd-even" ef- 
fect in the iron-peak element production. This shows [Mn/Fe] 
<~ -0.4 for [Fe/H] < 1 and a linear increase for [Fe/H] > 1 
up to [Mn/Fe] - 0.1 for [Fe/H] ~ 0.3 (Mc William, Rich, & 
Smecker-Hane 2003). 

Independent of the Mn II detection, the effects of dust ob- 
scuration can be estimated by assuming that this is propor- 
tional to the column density of Fe II. We know that in the so- 
lar neighborhood, the optical obscuration due to a gas cloud 
is (Bohlin, Savage & Drake, 1978): 



gives, for SFR = 1 M© yr luminosities between 8.7 and 7.2 
x 10 27 erg s" 1 Hz" 1 for the 1500-2800 A interval. We adopt 



(6) 



As ZNhi Npeii, Eq. 6 can be expressed in terms of the Fe II 
column density in the cloud: 



Ay ~ 0.5 x io losA,MI+5Fe " 21 " log ( Fe / H ) 



(7) 



where log(Fe/H) = -4.50 is the solar iron abundance rela- 
tive to hydrogen. For a Fe dust depletion 5p e — 1 dex (similar 
to the warm disk and warm disk+halo depletion in the MW; 
Savage & Sembach 1996), and logAWi = 15.54^, we ob- 
tain Av = 0.6+Q J, regardless of the H I column density. At the 
wavelength of the Ha the extinction is <~ 0.5 magnitudes, or 
a factor of <~ 1.6 in the flux. This is the value used in §4 to 
derive SFRs from [O II] emission. For <5p e = 0.6, like in the 
MW warm halo clouds, the optical extinction is Ay = 0.2+q 2 . 

In QSO-DLAs (were iron depletion is on average fe = 
0.52 ± 0.23) the optical extinction is generally lower 
(Ay ^0.1), most likely due to an observational bias (Savaglio 
2000). This is marginally supported by the work of Ellison 
et al. (2001), who tentatively found more DLAs in radio- 
selected QSOs than in optically selected QSOs (~ 2c signifi- 
cant). 

8.3. SFRs from UV luminosities 

Independently of the [O II] luminosity, SFRs can be derived 
using the rest-frame UV luminosities. These are more af- 
fected by dust extinction correction than those obtained us- 
ing [O II]. Moreover, UV emission is less directly connected 
to the instantaneous star formation, because associated with 
less massive stars than the nebular emission. However, it is 
interesting to evaluate SFRs from a completely independent 
method. TheM2ooo absolute magnitude (Table 1) gives an ab- 
solute intrinsic (dust corrected) luminosity: 

i 2 ooo = 4irR 2 Wpc x 10- 19 - 438 -°- 4 ( m ™o- a ™o) (8) 

where R\o pc = 3.08 x 10 19 cm, and A 2 ooo is the global (in the 
gas and stars) extinction at 2000 A. To derive SFRs from UV 
luminosities, different conversion factors can be used, accord- 
ing to different models that uses different initial mass func- 
tions (IMFs) and metallicities. Even if the conversion factor 
can vary from one model to another, Glazebrook et al. (1999) 
have shown that the variation is at most within a factor of two. 
For convenience, we select an intermediate case valid for a 
solar metallicity and Salpeter IMF. Glazebrook et al. (1999) 



SFRzoootMoyr" 1 ) 



i-2000 



10 27 < 



Hz" 



(9) 



To derive SFR2000, we nee d the dust extinction at 2000 
A (Eq. 8). We use two approaches: the MW extinction law 
and the Calzetti law (Calzetti 2001) derived from a sample of 
local starburts. The MW extinction is more appropriate for a 
uniform foreground screen of dusty gas (with covering factor 
1). In this case, using Ay = 0.6 obtained from the Fe II col- 
umn density, we derive Eg-v = 0.19. On the other hand, if the 
covering factor of the gas is < 1, then the dust extinction is 
patchy and differs on average from the one derived from the 
dust depletion in the gas. In this case we use Eg-v = 0.44, the 
mean value derived from a sample of local starbursts (Calzetti 
et al. 1994). For both cases we obtain A2000 ~ 1.7 magnitudes. 
SFRs are given in Table 1 ; the mean (and median) value is 36 
M yr -1 . For the subsample of galaxies with measured [O II] 
emission, the mean is 32 M yr -1 , close to the 40 M yr -1 
derived from the [O II] emission. A <~ 1 .7er significant linear 
correlation is found between the two sets of star formations 
(with SFR[ Q11 ] <~ 1 .2 SFR2ooo)- Our SFR2000 is measured from 
the aperture corrected photometry, that includes the UV flux 
from a large part of the galaxy. The consistency with values 
derived from the [O II] observed in the spectra (these are as- 
sociated with the smaller H II regions) suggests a high degree 
of mixing between the gas and the stars. 

8.4. Spectral slope in GDDS galaxies 

If the flux can be approximated by the power law f\ °< \@, 
we can estimate the spectral slope (3 in the individual spectra. 
The observed spectral range is not very large, 5000-9800 
A, therefore the spectral slope measured from there can be 
inaccurate, also because the flux calibration in MOS data is 
not perfect. However, we can use the observed Vega V and 

1 magnitudes to derive (3. If the AB magnitudes are Vab = 
V + 0.02 and I AB = 7+0.45, then (3 is: 



a /-V + 0.43 „ 

/3 = 0.4x— — — -2 

log(Ay/A/) 



(10) 



where Ay = 5500 A and A/ = 8000 A. For a galaxy at z = 1 .6, 
(3 would measure the spectral slope in the rest-frame inter- 
val 2100-3100 A. For the individual 13 GDDS galaxies we 
obtain a mean value ((3) = -1.15 with a large dispersion of 
0.9. In local star-forming galaxies, Calzetti et al. (1994) 
have measured (3 in a sample of 39 local starbursts in the 
wavelength range 1250-2600 A. From this sample we de- 
rive a mean value ((3) = -1.05, with a dispersion of 0.65. 
At high redshift Shapley et al. (2003) derived (3 in LBGs us- 
ing the intrinsic UV colors (A — 1200- 1800 A) and found 
-1.1 < (3 < -0.7. To compare our spectral slope to those 
found at lower wavelengths, we computed spectra using stel- 
lar synthesis models. We assumed a simple stellar population, 
with Salpeter IMF and continuous star formation. The spec- 
tral slope in the range A ~ 1200-2000 A is slightly flatter 
than in the range A <~ 2000-3000 A, due the contamination 
by strong stellar absorptions at A = 1200-2000 A. 

By applying the MW and starburst dust extinction laws, as 
described in the previous section, we calculated an intrinsic 
spectral slope in the interval A ~ 2 1 00 - 3 1 00 A for our GDDS 
galaxies (3i = -2.1 and -2.0, respectively. 
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Dunne, Eales, & Edmunds (2003) have found from submil- 
limetre surveys that dust content was much higher in the past, 
therefore most of the metal rich systems cannot be found in 
surveys that are too shallow, or too sensitive to dust obscura- 
tion. GDDS cannot detect galaxies or regions of galaxies that 
are strongly obscured, so it is certainly affected by dust, but 
probably much less than other galaxies surveys, because it is 
unbiased and highly complete for galaxies with K < 20.6. 

9. SUMMARY AND CONCLUSIONS 

The GDDS (Abraham et al. in preparation) is a public sur- 
vey that uses GMOS at Gemini, combined with the Nod & 
Shuffle technique, to study the properties of galaxies in the 
redshift desert (0.8 < z < 2.0). The GDDS observed galaxies 
with Iab < 24.8 and K < 20.6, and is the most sensitive survey 
focused on the redshift desert obtained so far. The K20 sur- 
vey obtained at the Very Large Telescope (VLT) is aiming at a 
similar scientific objective (Cimatti et al. 2002a) with 82 spec- 
troscopically confirmed galaxies in the redshift desert in 52 
arcmin 2 . In the first two GDDS fields observed, 60.5 arcmin 2 
and 59.5 hours exposure time, we identified 68 galaxies with 
spectroscopically-confirmed redshift. Among these, 34 are at 
1.3 < z < 2.5 (0.56 galaxies arcmin -2 ), whereas the K20 (17 
VLT nights) identified 17 galaxies (0.32 galaxies arcmin -2 ) 
in the same redshift interval, showing the high efficiency of 
GDDS. 

In this work, we have selected the UV brightest (AB abso- 
lute magnitude M2000 < -19.65) subsample with 13 galaxies 
that are in the redshift range 1.34 < z < 1.97 ((z) = 1.60) in 
order to study heavy elements in the neutral ISM. The selected 
galaxies are the most massive in the sample (M ^ 10 10 M©) 
and represents between 25 and 38% of the total number of 
galaxies found in the z > 1 .3 range (depending whether galax- 
ies with tentative or photometric-only redshifts are included 
or not). Only three of these 13 galaxies have K < 20.0. 

We detected strong [O II] emission associated with H II re- 
gions, with absolute luminosities ranging 5 - 40 x 10 41 erg 
s -1 , and roughly estimated SFR = 13-106 M© yr -1 and 
<SFR>~4OM yr -1 . The mean value is close to the mean 
value derived using UV luminosities found from the aperture- 
corrected photometry (32 M yr -1 ). This supports the as- 
sumption that on average the gas in the H II regions and the 
UV emitting stars are well mixed. 

The stack of the 1 3 spectra shows Fe II, Mg II, Mn II, and 
Mg I absorption. The column density of Fe II and Mn II 
provide a good measure of the total column density of iron 
and manganese in the neutral ISM. We found logA% e ii = 
15.54+jJ; 2 f and log/V M nii = 13.49 ±0.13. These values are 
larger than what typically measured in absorption systems 
along QSO sight lines (QSO-DLAs). Only 2 out of 24 
at 1.0 < z < 2.0 have similar Fe II column densities. On 
the other hand, high column densities are found in other 
star-forming objects. The Lyman break galaxy cB58 has 
log7V FeII = 15.25, while GRB host galaxies (GRB-DLAs) 
have logA^Fdi = 14.8-15.9. From the depth of the Mg II ab- 
sorption, we derived a conservative lower limit to the ISM 
filling factor of 70%. 

The high filling factor of the ISM and the consistency be- 
tween SFRs derived from [O II] and UV emission, suggests 
a homogeneous distribution of stars and gas. If we assume 
that the mean H I column density in the ISM is \ogNm < 
21.7 (as in 98% of QSO-DLAs), the mean metallicty would 
be Z/Z > 0.2 for a moderate Fe dust depletion correction 
(S Fe — 1 dex). High metallicities (Z/Z© = 0.5-2) are also 



found in H II regions of galaxies at 0.3 < z < 0.9 (Kobulnicky 
et al. 2003; Lilly et al. 2003). 

From the Fe II column density, we estimated an optical dust 
obscuration Ay ~ 0.6 mag. The mean spectral slope in our 
GDDS galaxies around the 2100-3100 A interval is (3 ~ 
-1.15, from which we derive an intrinsic dust corrected mean 
spectral slope /3 <~ -2. 

When completed, the GDDS will provide observations of 
two additional fields. From the full data set we will detect or 
give significant limits on other absorption lines. For instance, 
the Zn II detection at ~ 2000 A, or Si II at A ~ 1800 A, are 
powerful diagnostics for metal enrichment and dust depletion. 
GDDS galaxies and GRB-DLAs are clearly probing high-z 
galaxies with a strong star formation activity, for which metal 
enrichment and/or dust obscuration can be higher than in the 
QSO-DLA population. Our conclusion is that direct detec- 
tions of galaxies, and indirect through QSO absorption lines, 
provide a more complete picture of the galaxy census in the 
high redshift Universe. The former can be seen only if star 
formation is high, the latter only when metal enrichment is 
low. 
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